Probing Unified Origin of Dark Matter and Baryon Asymmetry at PAMEL A/Fermi 
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We propose an unified model of dark matter and baryon asymmetry in a lepto-philic world above 
the electroweak scale. We provide an example where the inflaton decay products subsequently 
generate a lepton asymmetry and a dark matter abundance with an unique coupling in the early 
universe, while the present day decay of the dark matter through the same coupling gives rise the 
observed cosmic ray anomalies at PAMELA and Fermi Large Area Telescope. 
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The recent observed anomalies p], [H, H[ in the cos- 
mic rays have given a lot of excitement. It has been 
shown that there is a clean excess of absolute positron 
flux in the cosmic rays at an energy E > 50 GeV 
even if the propagation uncertainty [5j in the secondary 
positron flux is added to the Galactic background. This 
leaves enough motivation for considering particle physics 
motivated dark matter (DM) models, such as annihila- 
tion @,0| or decay .7, § of DM, as the origin of positron 
excess in the cosmic rays 1 . However the origin of the 
DM, being interpreted as a long lived particle, goes be- 
yond the standard model (SM). Moreover, the origin of 
observed matter antimatter asymmetry and the origin of 
inflation are also two crucial phenomena which require 
the physics beyond the SM. 

In particular once cosmic inflation occurs, then af- 
ter inflation it must pave the way to excite not only 
the observed SM quanta but also the DM. This can be 
achieved minimally if the inflaton, <fi, itself carries the 
SM charges as in the case of the following examples 
which relies on the minimal supersymmetric SM (MSSM) 
setup 0, [H 03 • In the MSSM case the inflaton de- 
cays into the SM quarks and leptons and through ther- 
mal scatterings the lightest supersymmetric DM particles 
were created which are absolutely stable. On the other 
hand if inflation does not happen in the observable sec- 
tor, for instance if it belongs to a hidden sector, as in the 
case of plethora of examples [13| , then the onus will be to 
explain how to generate the desired degrees of freedom, 
i.e. SM baryons and dark matter abundance. 

The aim of this paper is to illustrate an example where 
lepton asymmetry [14| and DM abundance are generated 
right above the electroweak scale, i.e., 0(100) GeV. Our 
building block of beyond the SM physics is based on a 
crucial observation of cosmic ray anomalies 0, which 
may unravel the mystery of these issues in a unifying 
framework where the DM particle itself carries a net B—L 
asymmetry and decays very slowl y to the SM particles. 
Note that previous attempts [H, Tl6l H3 were made to 
unify dark matter and baryogenesis, but it is more chal- 
lenging to address why the DM annihilates/decays pri- 
marily into leptons and anti-leptons in a unifying frame- 



work, thus explaining the observed cosmic ray anomalies 
at PAMELA [l[ and Fermi @. 

For the purpose of illustration, let us augment the SM 
by adding a new £/(1)b-l gauge symmetry, and without 
supersymmetry. The anomaly free gauged B — L symme- 
try then naturally accommodates a new fermionic dark 
matter Al(1, 0,-1), where the quantum numbers inside 
the parenthesis shows the transformation properties of 
N L under the gauge group SU(2) L x U(1)y x C/(I)b-l- 
At a high scale the U(1)b-l gauge symmetry is bro- 
ken by a scalar field and give mass Mm = F«b-l to 
N = (1/V2)(N L + N C L ), where "v B -l" is the vacuum 
expectation value (vev) of the L t (1)b-l breaking scalar 
field which carries B — L charges by two units and F is 
the coupling between B — L breaking scalar field and Nl- 

In the broken phase of the £/(1)b-l gauge symme- 
try the Lagrangian involving the interactions of Nl, 
and the new massive charged scalars, rf~ (1,-2,0) and 
X~(l, —2,-2), can be separately written in terms of B—L 
conserving and B—L violating parts. The B—L violating 
part of the Lagrangian is given by: 



1 For astrophysical origins, see Ref. [j| and references therein. 



£a(b-l)#o = ^{M N ) a0 {N aL ) c Np L +m z rf x+h.c. , (1) 

with m 2 = //ub-Li while the relevant B — L conserving 
part of the Lagrangian is given by: 

+ K p rfN~L~e m + Upx^aLipL + h.c(2) 

where the indices a, (3 = e, [l,t represent the flavor basis 
of the SM fermion fields. In equation ©, £ L (2, -1, -1) 
and £r(1,— 2,— 1) represent the lepton doublet and sin- 
glet respectively, while H\(2, 1,0), #2(2,1,0) are two 
doublet Higgses which couple to up and down sector of 
SM fermions. 

An important point to note is that the mass term, 
m 2 rfx, which violates B — L by two units, gives rise 
a mixing between n and x- Due to this mixing the 
Majorana fermion N will be an unstable leptonic DM. 
Therefore, the decay products of N are only SM lep- 
tons/antileptons. If the lifetime of N is about O(10 25 ) 
s or so, then the decay products of N can naturally 
account for the observed excesses at PAMELA and 
Fermi. This observation will place non-trivial constraint, 
such as m « M V ,M X . Further we note that the only 
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coupling u h" is responsible for the production of a net 
lepton asymmetry and DM in the early universe, while 
the three body decay of DM through the same coupling 
at current epoch gives rise to the observed anomalies at 
PAMELA and Fermi. 
Baryon asymmetry: 

Let us now consider how can we explain the observed 
matter-anti-matter asymmetry in our setup. A natural 
possibility is to generate a lepton asymmetry [l?], HI] 
from the out-of-equilibrium decay of ?y~ field. Similar 
to Dirac leptogenesis [l|| , here we will argue that the re- 
quired baryon asymmetry can be generated from a con- 
served B — L number. In order to generate the baryon 
asymmetry we need following three steps: 

1. At first the CP-violating out-of-equilibrium decay 
of r\~ must generate an equal and opposite B — L 
asymmetry between Nl and £r. These two asym- 
metries should not equilibrate above the electro 
weak (EW) phase transition. 

2. Above the EW-phase transition the B — L asym- 
metry stored in £r gets transferred to £l, while 
keeping an equal and opposite B — L asymmetry in 

3. The B — L asymmetry stored in II then gets con- 
verted to a net baryon asymmetry in the presence of 
SU(2)l sphalerons, while keeping the B — L asym- 
metry stored in Nl intact. 

Note that all three steps happen right above the EW 
scale. Since r]~ is neutral under B — L, it can decay to 
a pair of lepton (£r) and antilepton (Nl)- Note that the 
decay of rj~ can not produce any lepton asymmetry since 
its decay does not violate any lepton number. However, 
if there are at least two rf~ fields, say rjf and rjf then 
there can be CP violation in the decay of rj~ fields. In 
their mass basis, spanned by ipf and ip^, the lightest ip 
field, say ipf , can generate a net CP asymmetry through 
the interference of tree level and self energy correction 
diagram [2(| ■ The CP asymmetry is then given by 
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Now assuming 



M„, M, 



-0(1), '''[-, 0(1) 



M m A1 V2 



and h}< * h% = O(10- 2 ) , (5) 



we get from Eqns. (J3J), Q and (J5J) the CP asymme- 
try ei ~ 10 -5 . Due to the CP violation the decay of 



-01 generates an equal and opposite B — L asymmetry 
between Nl and £r. Since the interaction between Nl 
and £r through the coupling 'h' is already gone out-of- 
equilibrium, the asymmetries between them don't equili- 
brate any more at the required scale of 0(100) GeV. On 
the other hand, the lepton number conserving process: 
£r£ c r £l£li mediated via the SM Higgs, remains in 
thermal equilibrium above the electroweak phase tran- 
sition. As a result the B — L asymmetry stored in £r 
gets transferred to £l through this L- number conserv- 
ing process, while leaving an equal and opposite B — L 
asymmetry in Nl- The transportation of B — L asym- 
metry from in to II can be understood as follows. Let 
us define the chemical potential associated with the £r 
field as jjL e R = Mo + Mbl, where /jbl is the chemical 
potential contributing to B — L asymmetry and (1q is 
independent of B — L. Hence at equilibrium we have 
the chemical potential associated with £l is given by 

MeL = f-eR. + = MBL + MO + M-ff- Tnus we see tnat 

the same chemical potential is associated with £l as of 
£r. Therefore, the net B — L asymmetry stored in £r 
can be passed on to £l- Since the SU(2)l sphalerons 
are in thermal equilibrium at a scale above 100 GeV, the 
B — L asymmetry stored in £l can be converted to a net 
baryon asymmetry, while an equal and opposite B — L 
asymmetry will remain in Nl ■ The two asymmetries will 
equilibrate when Nl will decay through the B — L vio- 
lating process. The net B — L asymmetry thus produced 
can be given as: 



3 T R 
Vb-l = - A B v ei 

4 rutf, 



(6) 



where B v is the inflaton branching ratio, which is of or- 
der, 0(1), and T R > 100 GeV is the reheat temperature 
of the universe and is the inflaton mass. 

The conversion of lepton asymmetry to the baryon 
asymmetry is obtained by tjb — {28/79)t]b-l- For 
Ti?/m w 10~ 4 and e x w 10~ 5 , we can achieve the ob- 
served baryon asymmetry t\b ~ 0(1O~ 10 ). A crucial 
point to note here is that the lepton asymmetry is virtu- 
ally created by a non-thermal decay of the inflaton decay 
products, r\ and Xi which we will discuss below. 

An obvious danger of washing out this asymmetry 
comes from the B — L violating process Nl£r — > £l£l 
through the mixing between rj and \- However, this 
process is suppressed by a factor (m 2 /M 2 M^) 2 for 
to <C M V ,M X , and hence it cannot compete with the 
Hubble expansion parameter at Tr ~ 100 GeV. Hence 
the net B — L asymmetry produced by the decay of 
77 will be converted to the required baryon asymmetry 
without suffering any washout. 

Dark matter abundance : 

Next we discuss the number density of the lightest Nl to 
check if it satisfies the observed DM abundance. It turns 
out that a large abundance of Nl will be produced non- 
thermally by the decay of 77, which is also non-thermally 
produced by the inflaton decay. 
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As we will argue below, the annihilation cross sec- 
tion of Nl is larger than the canonical one, (cr|u|) c ~ 
3 x 10 -26 cm 3 s~ 1 . In this case the final abundance of 
the thermal component is much smaller than the obser- 
vational value, 



Ydm = = 4 x 10 
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1 TeV\ ffl-ouh 2 



Mdm ) 



0.11 



(7) 



where Mdm is the DM mass, Hdm is the density pa- 
rameter of the DM with h being the normalized Hubble 
constant. 

Since we consider a case when the inflaton, 0, decays 
well after the standard freeze-out epoch of the thermally- 
produced Nl, i.e., at cosmic temperature T < Mat/25. 
Then the yield value of Nl is estimated by 



Yj 



n N , 



Nl 



3 T R 

4 'met, 



(8) 



For Tfi/m^ ~ 10 -4 , as required by the lepton asymme- 
try we find that the relic abundance of Nl is given by 
Ym l ~ 10~ 5 , which is much larger than the observed DM 
abundance (J7J). 

However thanks to the gauge coupling of Nl, such that 
Nl can now annihilate into the SM fermions through 
the exchange of Zb-l gauge boson after its non-thermal 
production. We can then obtain the final abundance of 
Nl by solving the Boltzmann equations: 
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(9) 




m/T D x=m/T 



FIG. 1: We illustrate the time evolution of pdm/s as a func- 
tion of x = Mdm/T, where pom/s = Mdm^dm- %f is the 
thermal- freezeout value of x (~ 25). To is the cosmic tem- 
perature when DM is non-thermally produced (To = Tr 
in the current model). The curved line traces the thermal- 
equilibrium value of DM when the thermal-production pro- 
cess is important, i.e., x < xf, which was omitted in Eq. 
The diagonal thick (thin) solid line represents freeze-out value 
when (<t|u|) is larger than (or equal to) the canonical value. 
The upper-horizontal dashed line shows the observational 
abundance. In the current model, the thermal component 
(lower-horizontal dashed line) is smaller than the observa- 
tional value. The star symbol represents the point when the 
initial value was put at T = To — Tr. Just after this non- 
thermal production, the DM annihilates immediately, and its 
abundance is reduced to the crossing point between the diago- 
nal thick-solid line and the vertical line, from then it becomes 
constant. 



where (cr\v\) as (1/Att)M^/v b _ l , and we have omit- 
ted the production term from the thermal bath, 



(a\v\)r 



'Nl ,eq 



in the right-hand side of the second line. 



Then we approximately obtain, 



Y, 



N L 



{(TV S 



(10) 



The right-hand side of Eq. (jTUJ) is oc 1/Tr, which means 
that the late-time decay induces a larger freeze-out value. 
In Fig. [1] we illustrate such a non-thermal production 
and/or further annihilation mechanism. 

By equating Yn l ~ Ydm we get a constraint on the 
B — L breaking scale with satisfying the observational 
DM density to be 



w B -l - 5 x 10 3 GeV 



DM' 



0.11 



1/4 



100 GeV 



M N \ 
3 TevJ 

1/4 



1/4 



(11) 



This gives (a\v\) - O(10~ 25 ) - O(10- 24 ) cm^" 1 for 
Mjv ~ O(TeV), which is larger than the canonical anni- 
hilation cross section and makes the thermal component 



sub dominant. In turn, we understand this mechanism 
intuitively by using Fig. [1] as follows. If we specify a 
U(1)b-l breaking scale (~ TeV in this model) or an 
annihilation cross section, the model necessarily has a 
crossing point between the diagonal and the horizontal 
lines at x = Mdm /Tr shown in Fig. [TJ which gives the 
right observational abundance of DM and its production 
epoch. In the current model, we demand the cross sec- 
tion to be larger than the canonical value and Tr to be 
larger than 100 GeV. 

Note that the annihilation process NlNl — > // is a 
B — L number conserving process and therefore does not 
transfer any B — L asymmetry to the SM fermions. As 
a result the B — L asymmetry produced via the decay 
of will survive until far below the electro-weak phase 
transition. 

Since B — L is already broken, the lightest Nl is no 
more stable. It will decay through the three body pro- 
cess: Nl — > e aR e ^L^jL, with /3 7, through the mixing 
of r\ and \- Since the coupling of x to two lepton dou- 
blets is antisymmetric, i.e., (3 ^ j, the decay of Nl is not 
necessarily to be flavour conserving. In particular the 
decay mode: N L -> r^r^T7 eL (l7 flL ) ] violates L e (L M ) by 
one unit while it violates L = L e + L„ + L T by two units. 
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In the mass basis of Nl the life time can be estimated to 
be 



0.1 



tat = 2.0 x l(T J s 



50 GeV 



26- A° 



-2 X 2 / 10 -7N 2 
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m \ s / 1 TeV 



10 6 GeV J V m n 



(12) 



where we assume that M v ~ M x « in order to get 
a lower limit on the lifetime of Nl- The prolonged life 
time of Nl, may explain the current cosmic ray anoma- 
lies observed by PAMELA [l[ and Fermi @|. The elec- 
tron and positron energy spectrum can be estimated by 
using the same set-up as in Ref. ||. In Figs. [2] and 
|3] we have shown the integrated fluxes in a typical 
decay mode: Nl ^ t~t + D up to the maximum avail- 
able energy Mjv/2 for tn = 4.0 x 10 25 sees. From there 
it can be seen that the decay of Nl can nicely explain 
the observed positron excess at PAMELA and e^- ex- 
cesses at Fermi. While doing so we assume that the 
branching fraction in the decay of N — L to t~t + v is 
significantly larger than the other viable decay modes: 
Nl — > n~ n + v and Nl — > e~e + v. However, if the decay 
rate of Nl — > is comparable to Nl — > t~t + D 

then it can explain the observed anomalies at PAMELA 
and Fermi, while the decay mode: Nl 



i e + v produces 
fluxes at Fermi and therefore unfavorable. 
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FIG. 2: Positron excess from Nl — > t t v with Mn = 3 
TeV. The fragmentation function has been calculated using 
PYTHIA [H. 



Another potential signature of this scenario is the 
emission of energetic neutrinos from the Galactic cen- 
ter [25| which can be checked by future experiments 
such as IceCube DeepCore [H and KM3NeT_j27|. We 
will come back to this issue in a future work 2£ 



> 
O 

+ 

w 



0.01 



0.001 



1-crmi 
HESS:08 
HESS:09 
ATIC: 

PPP-BETS:0S 




Background 



10 



100 

E(GeV) 



1000 



FIG. 3: e excess from Nl — * T T V with Mn = 3 
TeV. The fragmentation function has been calculated using 
PYTHIA [3. 



Inflation and reheating: 

So far we have not discussed anything about inflation sec- 
tor. Let us now consider a hidden sector inflaton <fi(l, 0, 0) 
which excites the observed DM abundance and SM lep- 
tons during the process of reheating which occurs after 
the end of inflation. Let us consider a simple toy model 
of inflation, where the inflaton potential admits a point 
of inflection, see [10, fill fl5j| : 



v(4>) - > 2 



An 



+ 12< 



(13) 



where A w Am^ and k ~ 10 10 . Inflation can happen 
near a point of inflection given by </> ~ VSttl^/k ~ 
10 16 GeV with an Hubble expansion rate, Hi n f ~ 
(ml/nMp) - 10 4 GeV. The amplitude of the density 

perturbations will be give by: 6h ~ {l/5n)(Hf n} /<j)) ~ 
5 , where the number of e-foldings 
j22j |. One of the dynamical prop- 



(K 2 M P /3m,p)N 2 - 10 
is given by: N 2 ~ 10 3 
erties of an inflection point inflation is that the spec- 
tral tilt can be matched in a desired observable range: 
0.92 < n s < 1.0 for the above parameters [13, LH LH, El . 

The inflaton decays into heavy charged scalars 77 and 
X- Let B v be the branching fraction in the decay of <j> to 
x be the branching fraction in the decay of 4> to 
As we discussed above the charged scalars rf 1 and 



Tj^~ and B 



X 

X ± couple to the SM degrees of freedom. The reheating 
occurs when the inflaton begins oscillations. The largest 
decay rate happens for the largest amplitude of oscilla- 
tions, for instance when (<p) ~ <po, see for instance [23|. 
As a result the Universe gets reheated up to a desired 
temperature: 



T R - O.l/tVWpi = 1.2xl0 2 GeV 



(lO- 17 ) 



( m <f> \ 



1/2 

Vl0 6 GeVy 

(14) 
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above the electro-weak scale to facilitate a successful 
baryogenesis, where = (g 2 /8tt) (((/)} /m^) 2 is the 
decay rate of <j> with g is the quartic coupling: g<f> 2 \ry rj + 
X^X)- A small decay rate of inflaton to rj^ and x^ ensures 
the optimal temperature just right above the electro- 
weak phase transition. 

To summarize, we have explored a simple model of 
lepto-philic universe where DM carries a net B — L 
asymmetry and decays only into the SM leptons can 
explain the observed positron excess at PAMELA and 
excesses at Fermi. These anomalous observations at 
PAMELA and Fermi may indirectly probe the common 
origin of mysterious DM and baryon asymmetry as we 
have shown. The baryon asymmetry in our model is cre- 



ated via lepton conserving leptogenesis mechanism which 
gets converted into baryon asymmetry via the electro- 
weak sphalerons. Before closing we note that the model 
explained here can be embedded within supersymmctry 
without further challenges, neither the leptogenesis nor 
the DM mechanisms will alter, the parameters and the 
observed values of the inflationary perturbations and the 
tilt in the spectrum would remain so. The details will be 
published elsewhere [28j |. 
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